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Truncated Overlap Fermions 
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In this talk I propose a new computational scheme with overlap fermions and a fast algorithm to invert the 
corresponding Dirac operator. 



1. INTRODUCTION 

After many years of research in lattice QCD, 
it was possible to formulate QCD with chiral 
fermions on the lattice |l| Qj. 

The basic idea is an expanded flavor space 
which may be seen as an extra dimension with 
left and right handed fermions defined in the two 
opposite boundaries or walls. 

Let N be the size of the extra dimension, 
Dyy the Wilson-Dirac operator, and m the bare 
fermion mass. Then, the theory with Domain 
Wall Fermions is defined by the action 



Sdw ■= ^M^f 

+P+ip i+ i + P-tpi-i], 
P+Ojv+i + mipi) = 0, 



mtp N ) 
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where M. is the five-dimensional fermion matrix 
of the regularized theory and Z?l = M — Dw with 
M £ (0, 2) being a mass parameter. 

In this talk I define a theory with Truncated 
Overlap Fermions in complete analogy with the 
Domain Wall Fermions by substituting 



P+^i+x 
P-^-i 



(£>H + l)P+^+i 
(£>"+ l)P-^i- X 



(2) 



while the boundary conditions remain the same 
as before. 

Both theories can be compactified in the walls 
of the extra fifth dimension as low energy effective 
theories (see below) with the chiral Dirac opera- 
tor D satisfying the Ginsparg- Wilson relation || : 

jsD- 1 + prHs = 2aj 5 R, (3) 

where a is the lattice spacing and R is a local 
operator trivial in the Dirac space (see below for 
R- locality tests). From now on I set a = 1. 



I defined Truncated Overlap Fermions such 
that in the large N limit one obtains Overlap 
Fermions j| with the Dirac operator given by Q : 



Dov = 



1 



1 — m 



7 5 sgn(if) 



(4) 
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where H = 75-D" . 

Until now computations with chiral fermions 
and standard algorithms have been very expen- 
sive. The extra fermion flavors introduce a large 
overhead. One multiplication with the fermion 
matrix costs O(n) Diy-multiplications with n ~ 
N for Domain Wall Fermions and much larger for 
the overlap operator . 

In this talk I propose a fast algorithm which 
makes these simulations an order of magnitude 
faster. The key observation is the lack of gauge 
connections along the fifth dimension. 



2. TRUNCATED 
FERMIONS 



OVERLAP 



I recall the action of the Truncated Overlap 
Fermions: 

Stov :=M-D" -1)^i 
+(£)!! + l)P+tp2 - m(L>H 

+(£>H + l)P+Vi+i + {D^ 
+$ N [{D\\ - 1)Vjv 

-m{DW + l)P+ipi + (£>N + 1)P_Vjv-i] 

Let P T be the matrix representing the unitary 
transformation: 



1)P-Vjv] 
1)P-^-i] 



(5) 



(6) 



Xl = P+1pl + P-I^N, 

Xi=P+ipi+P-ipi-i, i = 2,...,N 

and S the matrix representing the diagonal trans- 
formation: Xi = i>il5{H — 1), i = l,...,N. Let 



Figure 1. Norm of D kernel in spin and color 
space with the distance r from the origin for N = 
4 (circles) and N = 32 (crosses). 



Figure 2. Norm of R kernel in spin and color 
space with the distance r from the origin for N = 
4 (circles), N = 32 (crosses) and N = 64 (pluses). 



also the transfer matrix along the fifth dimension 
be defined by: T = ±±§ . 

In the new basis I obtain the following action: 



(7) 



Stov = Xi[{P+ - mP_)xi - T X2] 
+ J2 i =i 1 Xi{Xi -Txi+i) 
+XN[XN-T{P_-mP+)xi] 

Integrating over the Grassmann fields I get: 

AetM = det[(P+ - mP_) - T N (P- - mP+)\ (8) 

where I ignore the Jacobian factor coming from 
the diagonal transformation. 

If All is the same matrix as M but with the 
special choice m = 1, I define the effective low 
energy theory with the Dirac operator given by 
the equations: 

D = {P T Mi 1 MP) 1A , D- 1 = (P T M- 1 M 1 P) 1 . 1 (9) 

where the subscript 1, 1 stands for the (1,1) block 
of an NxN partitioned matrix along the fifth di- 
mension. 

In terms of the transfer matrix the Dirac oper- 
ator can be written as: 



_ 1 + TO 1 - TO 1 — T N 

2 + "I -75 1 + T N 



(10) 



I can repeat this derivation for the Domain 
Wall Fermions with the sole changes Xi = 
$a h (HP + -\), i=l,...,JVandT= I= £ J!q: (l+ 
HP-), the rest of the formulae remaining the 
same. 



3. LOCALITY AND GINSPARG- 
SON RELATION TESTS 



WIL- 



I test the Ginsparg- Wilson relation for the 
Truncated Overlap Fermions, so that they can be 
used like Domain Wall Fermions. 

I computed the norm of D and R kernels in 
spin and color space with the distance r — ^Jx^x^ 
from the origin on 4 4 lattices at j3 = 6.0. In Figs. 
1 and 2 I show the maximum values at on-axis 
distances. Note that x M , \i = 1, 2, 3, 4 is measured 
modulo lattice size in that direction. 

Fig. 1 suggests exponential fall-off of ||Z?(r)|| 2 , 
whereas Fig. 2 shows that R approaches a 
Kronecker-Delta function as N grows. 

4. A FAST INVERSION ALGORITHM 

I use Truncated Overlap fermions to define the 
following 

ALGORITHMl (Generic) for solving the sys- 
tem Dqvx = b: 



Given N, x , r (= b), tol, toll, set tolo = 1 
and iterate: 
for i = l,... 
tol = tol toli 

Solve Dy = r ; _ x within tol 

Xi = + y 
r ; = b DovXi 

j/lkilb < tol, end for 



(11) 



where by o is denoted a vector with zero entries 
and toh, tol are tolerances, toh is typically orders 
of magnitude larger than tol such that the work 



per Dqv inversion is minimized. 

Remark 1. Bold face equations represent the 
smaller system soultion and the correction of the 
right-hand side. The straightforward application 
of the ALGORITHM 1 gives a two-level algo- 
rithm. By calling it again in solving the smaller 
system and iterating, one gets a multi-level algo- 
rithm. 

In Fig. 3 I compare the norm of the residual 
n = b — DovXi of the Conjugate Residual (CR) 
algorithm (which is optimal since Dov is normal 
(lO)) and ALGORITHM 1 . I gain about an order 
of magnitude (in average) on 30 4 4 -configurations 
at j3 = 6.0 and m = 0.1. For the coarse lattice I 
used N — 6 with the Truncated Overlap Fermions 
and the Lanczos method to compute Dov @j- 

Remark 2. Dynamical fermions can be im- 
plemented similarly. The corresponding Hybrid 
Monte Carlo (HMC) algorithm can be obtained 
by working with an approximate Hamiltonian in 
the coarse lattice and by a global correction on 
the fine lattice. 

One may use also as a starting point Truncated 
Overlap Fermions: detD =detA4/det.Mi. This 
way, all known simulation algorithms for dynam- 
ical fermions apply. 

5. CONCLUSIONS 

I showed that Truncated Overlap Fermions 
may be used in two ways: 

a) to implement Overlap Fermions in the same 
fashion as Domain Wall Fermions; 

b) to construct a multi-level inversion algo- 
rithm for Overlap Fermions which saves an order 
of magnitude of computer time compared to the 
state of the art methods. 

Further tests are needed to verify these results 
on larger lattices. 

Recently, the possibility of a Multigrid algo- 
rithm along all dimensions is raised |11|. In this 
case a gauge fixing is needed. 
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Figure 3. Norm of the residual error vs. the 
number of D\y multiplications on 30 configura- 
tions. Circles stand for the straightforward inver- 
sion and stars for the new algorithm. 
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